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Abstract 

The photoionization models for the narrow emission line regions of powerful radio galaxies have 
yet to be tested in depth. To this end, we present high-quality long-slit spectroscopy of the powerful 
double-nucleus radio galaxy 3C 321. The data have good enough spatial resolution to be able to trace 
the variation in emission-line properties on kpc scales. Continuum modelling and subtraction enables 
QO ' the faint emission line fluxes to be measured in several regions across the emission line nebula. We 

£\J , plot diagnostic line-ratio diagrams and compare them with the predictions of various photoionization 

£NJ 1 models, finding that the data is best fit by models which assume a mixture of optically thin and thick 

' clouds illuminated by a power-law continuum. 

The emission line kinematics, line ratios and deduced physical conditions show remarkably little 
variation across the source. We deduce a mean electron density of 400 ± 120 cm -3 and a mean 
temperature of 11500 ± 1500 K. Assuming a single population of optically thick line-emitting clouds, 
we calculate a mean photoionization parameter of (1.1 ±0.5) x 10~ 2 and hence a photoionizing photon 
<*S !' luminosity of Q ~ 10 55 — 10 56 phots _1 sr _1 . This indicates a central engine as luminous as that of the 

powerful quasar 3C 273, yet there is no evidence for such an energetically prolific central engine at 
either far-infrared or radio wavelengths. We therefore conclude that the mixed-media models, which 
give Q~ 5 x 10 53 — 5 x 10 54 , represent a more likely scenario. 

As a by-product of the continuum subtraction we infer that young stellar populations account for 
~ 0.4% of the visible stellar mass in the galaxy, and that these populations are spatially extended. 
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1 Introduction 

Emission lines are the most easily measured features of radio galaxies, and yet the physical mechanisms 
involved in their production are not well understood. Various suggestions have been made as to the 
method of ionization of the line-emittin g gas, including shocks from interactions between the radio jet s 
and the warm clouds ( Clark et al. 199§| ) and photoionization by a central source ( Robinson et al. 1994 ). 



Central source photoionization is particularly attractive in light of unified schemes for powerful radio 
sources, which hold that radio galaxies and quasars are a single class of object viewed from different 



directions ( [Barthel 1989| ; [Lawrence 199l| ). If these orientation-dependent unified schemes are true, and 
radio galaxies contain quasar nuclei obscured from direct observation by a dusty torus, we expect to 
see the effects of the powerful AGN illuminating the surrounding regions. Assuming that central source 
photoionization is the dominant mechanism, the emission line spectra may be used to infer the properties 



of the central engine, such as the shape of the photoionizing continuum and its power (Tadhunter et al. 



1994; Morganti et al. 1991). Indeed, emission-line luminosity has been shown to correlate broadly with 



178 MHz radio power (Tadhunter et al. 1998; Rawlings & Saunders 1991), suggesting a link between the 



AGN continuum emission and the mechanism of radio-jet production. 
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One way of investigating th e properties of the central engine is by m odelling the emission line 
spectrum of the photoionized gas ( Robinson et al. 1987 ; Binette et al. 198g| ). In general, these models 
depend on both the properties of the line-emitting medium and the ionizing spectrum, which in turn 
depends on the properties of the central AGN. The traditional approach has been to assume a power-law 
or hot black-body continuum ionizing optically thick, solar abundance clouds of gas. This approach has 
been successful up to a point, but there still remain several outstanding difficulties in explaining the 
observed spectra of radio galaxies in terms of pure photoionization: 

• Discrepant line ratios. Some line ratios are problematic in that they are not c orrectly predicted 
by ph o toionization models in some objects. Th ese include [OHT|A50 07/A4363 ([Tadhunter et al. 



1989a|; pVIorganti et al. 1997b| ), HeIIA4686/H/3 ( [Binette et al. 1996| ) and [NII]A6583/[OII]A3727 
( Tadhunter et al. 1994 ) . At least some of these discrepancies can be explained in terms of departures 
from the assumptions of the simplest models. For example, the discrepant [OIII]A5007/A4363, 
HeIIA4686/H/3 ratios c an be modelled if th e photoionized medium is a mixture of optically thin 
and thick components ( Binette et al. 1996), while th e [NII]A6583/[OII]A3727 can be explained in 
terms of non-solar abundances ( Viegas fc Prieto 1992 ). However, the mixed- medium model requires 
specific column depths in the optically thin and thick clouds and cannot explain some of the more 
extreme line ratios measured in some sources. 

• Spatial variations in line ratios. These are often difficult to explain in terms of a central 
ionizing source. For example, the [NII]A6584/Ha and [SII]A6726/Ha line ratios in Cygnus A show 
no variation in value across both core and extended regions, despite the fact that [OIII] line ratios 
vary dramatically across the same regions (Tadhunter et al. 1994). 



• The effects of shocks. Line broadening and differences in the widths of lines of different ionization 
suggest shocks are also important, even in the n arrow line region (NLR) where there is little direct 
morpholog ical evidence for jet / cloud interaction ( Clark et al. 199S ; Villar-Martin et al. 199£ ; Capettj 
et al. 1999[ ). 



To test the photoionization models, high quality spatially resolved spectra of a predominantly photoion- 
ized radio galaxy are required. We have obtained such data for the powerful radio galaxy 3C 321. The 
advantages of this dataset are that it has wide spectral coverage, accurate flux calibration, low airmass, 
subarcsecond seeing, S/N ratios high enough to measure electron density in several spatial bins (a feature 
not shared by any other radio galaxy with published spectra), and a high enough S/N in the continuum 
to allow accurate modelling and subtraction. 

The double nucleus host galaxy of 3C 321 is an ideal object for an in-depth study of photoionization 
models. Its relatively low redshift (z = 0.0961) means that ground-based long-slit spectroscopy allows 
the study of the variation of its emission-line spectrum on kpc scales. 20cm VLA radio images in 
combination with emission line images ( Baum et al. 1988| ) show little sign of jet-cloud interaction, despite 
the overall complexity of the emission line structure, and there is little evidence for the extreme emission- 



line kinem atics which would indicate the action of shocks, as observed, for example, in 3C 171 (Clark 



et al. 1998| ) . Polarization studies also provide evidence of illu mination by a cen t ral source, since scattered 
quasar feat ures are detected in both polarized and direct light ( Filippenko 1987 ; Draper et al. 1993 ; Young 
et al. 1996 ). Furthermore, the fact that 3 C 321 is at the brighter end of the 60 [im luminosity range 
of Narrow Line Radio Galaxies (NLRGs) ( Hes et al. 1995 ; Heisler fc Vader 1994 ) also suggests that its 
central region s contain a powerful illu minating source. In addition to the quasar component, continuum 
modelling by Tadhunter et al. (1996) (TDS96) showed evidence for a significant starburst component. 

We present the long-slit data for 3C 321 in Sections || & || and discuss the implications of these 
data for the photoionization models of powerful radio galaxies in Section |j. 



2 Observations 



2.1 Data acquisition and reduction 

The observations were made in July 1996 using the ISIS spectrograph on the 4.2m William Herschel 
Telescope on La Palma. Two sets of 1500s exposures were taken simultaneously on the red and blue arms 
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of ISIS giving a wavelength coverage of 3630-9820 A. Not all of this range was useful, as the dichroic in 
ISIS attenuates the region above ~ 5900 A in the blue spectrum, and flux calibration errors dominate 
above ~ 8200 A. Gaussian fits to the night sky lines show that the mean spectral resolution of our data 
is 8.2 ± 0.1 A in the blue spectra and 7.9 ± 0.9 A in the red. 

The seeing was measured to be 0.8 arcsec (FWHM), and the low airmass (1.008) ensured that the 
effects of differential atmospheric refraction were negligible — an important consideration given the large 
wavelength range covered by the data. The spatial scale of our spectra is 0.79 kpc per pixelQ. 

The slit was placed at a position angle of 130° on the sky, the direction of the axis defined by the two 
nuclei of the object. Figure [j] shows a diagram of the slit superimposed o n an [OHIJAA4959 + 5007 HST 
image taken with WFPC2 using a linear ramp filter ( Martel et al. 1999|) . A clear dust lane is visible in 
the SE blob, which is not resolved in our long slit images, although there is evidence for a significant 
increase in the E(B-V) reddening in this region (Region F in Figure - see Table |). Also marked are 
the 7 regions chosen for the detailed analysis described below. 

The data were reduced using the Starlink FIGARO package and IRAF. Bias subtraction was carried 
out and cosmic rays removed. The spectra were then wavelength-calibrated, corrected for atmospheric 
extinction and flux calibrated, using the standard stars and arc exposures taken at the same time as the 
science images. The data were also corrected for both S-distortion and tilt, which may be caused by the 
slight misalignment of the slit with the grating and CCD, by using the spectra of standard stars taken at 
the same time as the observations. After correction, the residual offset in the spatial direction between 
the two ends of the spectra was estimated to be < 0.3 pixels. The two blue and two red spectra were 
co-added and the composite blue spectrum was then re-sampled to give the same spatial pixel scale as 
the red, and the region of the spectrum severely attenuated by the dichroic was removed. Finally, the 
two composite spectra were aligned and shifted back to the rest-frame using the mean redshift measured 
from the brighter emission lines. 

The reduction process resulted in two spectra with useful rest wavelengths in the ranges 3314-5400 
& 5891-8200 A. We estimate a flux calibration error of 5% across this range. 



3 Results 

Some analysis was made of the data prior to the continuum subtraction. In particular, kinematics of the 
[OIIIJA5007 line were m easu red across the object on a pixel-by-pixel basis. The results of this line-fitting 



are presented in section |3.2 . 

Next, seven spatial bins were chosen for detailed analysis (see Figure |l|). These were selected to 
have some morphological meaning and/or to give a high enough signal-to-noise in the density diagnostic 
[SIIJAA6717, 6731 lines to measure densities 3<r away from the low-density limit. The HST image was used 
to identify spatial bins which fall either side of the central dust-lane, and to sub-divide the emission line 
and continuum spatial peaks whilst maintaining a high signal-to-noise (see Figure ||) . Note that Gaussian 
fits to the mean spatial profiles of the [OIII]A5007 and line-free continuum (A8849) show significant offsets 
between the continuum and emission line peaks in both nuclei, with the north-west emission line peak 
offset by 0.11 ± 0.04" (0.26 kpc) to the north-west of the north-west continuum peak and the south-east 
emission line peak offset by 0.66 ± 0.02" (1.6 kpc) to the SE of the SE continuum peak (Figure ^). 

3.1 Continuum modelling 

To enable the fainter emission lines to be measured more accurately, model continua were calculated 
and subtracted from the observed spectra. A nebular continuum was generated for each bin, comprising 
the blended higher (>H8) Balmer series together with a theoretical nebular recombination continuum 
(generated using dipso), both of which were normalized to the flux in H/3. T his was then convol ved with 



the H/3 line width and subtracted from the data prior to the modelling (see Dickson et al. 1995 ). 

Initially, the continuum in regions A-G was modelled using 2 components chosen to represent 
AGN light (either direct or scattered) and the light of the host elliptical galaxy: a power law of the form 
fx oc \ +a (corresponding to a frequency power- law, z/-( Q + 2 )) and a 15 Gyr old elliptical galaxy template 

1 Values of go = 0.5 and Ho = 50 km s -1 Mpc _1 are assumed throughout. 
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taken from Bruzual fc Chariot (1993)| . A more sophisticated model was then introduced, following TDS96, 
this time with an extra component corresponding to a contribution from a 0.1 Gyr starburst. The flux in 
36 continuum bins, selected from both the red and blue spectra, was measured for the data and the model 
components, avoiding emission lines, atmospheric absorption bands and image defects. The models were 
then generated by choosing a normalizing continuum bin and scaling the different model components so 
that the total model flux in the normalizing bin was less than 125% of the observed flux. In this way, a 
series of models was created, filling the 3- and 4-d parameter space defined by the relative fractions of each 
of the components, and the value of a. The best-fitting model was determined using a reduced %-squared 
test, with the fitting being performed by IDL routines written for this purpose. The contributions of the 
different components were varied between - 125% of the flux in the normalizing bin in increments of 
2.5%, and the power law was varied over the range —4.9 < a < 5.1 in 50 increments. The results of the 
modelling are shown in Tables @&|. These were used to generate and then subtract the best-fit model 
continuum from the observational data. The best fitting models for the 3-component fits are plotted in 
Figure together with the nebular-component-subtracted continua of the different spatial regions. 

We exclude region A from the discussion which follows. As Tables |j] & || show, the value of \ 2 
is much higher for this region than for the other regions. This is due the fact that the continuum in 
region A is both weak and noisy, and the background subtraction more uncertain. We therefore refrain 
from drawing any conclusions regarding the contributions of individual continuum components in this 
region, although we do subtract the best-fit continuum from the spectrum for the purposes of studying 
the emission-line spectrum. 

The 2-component fits are entirely adequate for the regions containing the two continuum peaks 
(C and F). For all other regions, however, the 2-component fits are poor, with the data exhibiting an 
excess of 10 — 15 percent of the continuum flux relative to the best-fit model in the wavelength region 
3750 - 4300 A. A similar discrepancy was noted by TDS96 for the integrated light of the nuclear regions. 
Moreover, the 2-component models also give an over-prediction of the flux above ~ 7000 A. Note that 
the increased red-coverage in our spectra compared with those of TDS96 gives us a 'longer lever' in terms 
of constraining the fits to the data. 

The 3-component models fit the data substantially better, giving much lower values of Xmin reduced 
and also giving a more convincing overall spectral shape. As suggested by the 2-component models there 
is no strong evidence for a young stellar population in regions C & F, but the other four regions have 
>20% of their continuum at 5060 A contributed by a starburst component. Note that the negative slope 
of the continuum between 6500 A and 9000 A, as measured in regions B and E (see Figure 3), is exactly 
as predicted for a significant young stellar population but is not expected for an old (15 Gyr) stellar 
population. 

These results are consistent with TDS96, whose 1-d data did not allow spatial resolution of the 
different regions presented here. Although they quote the contributions of the various components as a 
fraction of a bin centered at 3620 A, it is possible to estimate the contribution this would give at the 
5060 A bin we used. The TDS96 data yield contributions of 71, 21 and 8%, at 5060 A, of the continuum 
remaining after subtraction of the nebular component for the 15 Gyr, 1 Gyr and power-law components 
respectively, using their value of a = —1.7. (Note that their young stellar population is older than the one 
used in this paper.) Taking a broad spatial aperture corresponding as closely as possible to the one used 
by TDS, containing regions D - G, the continuum fitting gives contributions of 60, 20, and 17.5 % for the 
15 Gyr, 0.1 Gyr and power-law, which is entirely consistent given that the young stellar components in 
the two sets of models used have different ages, and also that our data provide much better constraints 
in the red. 

It is possible to calculate the total mass of stars present in each component from these results. 
The original [Bruzual fe Chariot (1993)| models are all normalized to a total mass of 1 Mq , so it is easy 
to scale the fluxes so as to give the contribution to the normalizing bin indicated by the best fit results 
in Table g. The masses thus derived are shown in Table ||. The young stellar population forms ~ 1% 
of the total stellar mass in all the spatial regions, and we believe this is the first published confirmation 
that the young stellar populations in radio galaxies are spatially extended. Note that for the older, 1 Gyr 
population used by TDS, we would expect the young stellar population to make a greater proportional 
contribution to the total mass. 

Although we believe our models show strong evidence for a young stellar population that is spatially 
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extended across the entire nebulosity, we emphasize that the models are not unique in terms of the age 
of the young stellar population. For example, TDS96 found that a 1 Gyr starburst provided a good fit 
to the data. More work is required to determine the age of the starburst accurately, although we stress 
that the emission line flux measurements are not significantly affected by the exact age of the starburst 



used in the continuum subtraction in the age range 0.01 — 1 Gyr (Dickson 1997). 



3.2 Emission line kinematics 

Plots showing the emission-line kinematics are presented in Figure |[ The bright [OIII]A5007 line was 
measured in every spatial pixel, whereas the other lines were measured in the regions identified above. 
Both velocity shifts and velocity widths of the lines are plotted, with the widths quadratically corrected 
for instrumental broadening using the measured widths of night sky lines as close as possible in wavelength 
to the emission lines being studied. The 'zero' of velocity shift is defined as the wavelength of the line 
measured at the centre of the continuum emission of the brighter (SE) nucleus, as determined by fitting 
a Gaussian to its spatial profile. 

The plots are encouraging in terms of using 3C 321 to test the photoionization models, since they 
suggest that, relative to other radio galaxies, 3C 321 is kinematically undisturbed. Not only are the 
lines relatively narrow (FWHM <^ 500 kms -1 ) across much of the nebula, but there are no significant 
differences in kinematics between lines of different ionization, as has been found in cases of jet-cloud 
interactions ( Villar-Martin et al. 1999j ). It is only in the farthest spatial regions to the SE and NW 



that there is any sign of broadening in the emission lines. In particular, region A (7.5 arcsec NW of the 
nucleus) shows evidence for significantly broader line profiles in all the measured lines (FWHM ~ 500—600 
kms -1 ). 

The spatial variation of radial velocities across the nebula is also small, with a spread of < 200kms~ 1 
in the range —6 to +2 arcsec. The variation which is observed may be due as much to the structure of 
the emission line nebula within the slit as to intrinsic velocity differences between different regions, since 
the seeing (~ 0.8 arcsec) was significantly less than the slit width. 

Generally, the values presented here correspond to typical narrow-line widths in radio galaxies 



1989a) 



(Baum et al. 199C), and are entirely compatible with purely gravitational motions (Tadhunter et al. 



3.3 Reddening 



The lin e fluxes were corrected for reddening using the peaton (1979) 



1975)| ex tinction curve, with Case B recombination line values from 



re-normalization of the [Nandy et al. 



Gaskell fc Ferland (1984)| for Ha/H/3 



and from Osterbrock (1989) for all other line ratios. The continuum-subtracted spectra allowed good fits 
to be made to H7 and US, as well as the other brighter members of the Balmer series, in all seven regions. 
The resulting E(B-V)values are shown in Table |^. These values suggest no intrinsic reddening in regions 
B & D, and a maximum in region F. The high values in regions E G are consistent with the HST image 
(Figure 0), which shows a prominent dust lane across the SE edge of region F. These values of E(B-V) 
are also comparable with the extinction measured in the dust lane of Cygnus A of E(B-V) ~ 0.3 — 0.6 



( [Tadhunter et al. 199J) . 

The estimated extinction values are smaller than the E(B-V) estimated by TDS96, but it should be 
noted that their measurement of Hy/H/3 was made before any continuum subtraction, and that therefore 
no correction was made for underlying Balmer line absorption. Their measured value of H7/H/3 was 
therefore ~ 25% lower than ours, which accounts for the apparent discrepancy in E(B-V). 

Despite their much larger value of E(B-V), TDS96 found that de-reddening the continuum had a 
negligible effect on the continuum fitting, in terms of the contributions of the different components. This 
suggests there is little to be gained in repeating the continuum fitting/subtraction for our spectra after 
correcting them for reddening. However, the emission line fluxes were themselves corrected for intrinsic 
reddening, using the values shown in Table [|. 
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3.4 Emission-line properties 

After subtracting the best-fit model continuum from the observed spectrum in each region, the emission 
line fluxes in each spatial region were measured using the Starlink package DIPSO and corrected for the 
effects of intrinsic reddening. The corrected fluxes relative to H/3 are shown i n Table p|. 



We can confirm the detection of a broad component to Ha reported by [Filippenko (1987)| in direct 
light. Whereas Filippenko only detected the broad Ha in the brighter nucleus (corresponding to apertures 
E & F), we also detect it in the fainter nucleus (aperture C), albeit at low signal-to-noise. The FWHM 
velocity- widths of this component, where measured, a re shown in Table |j. These are consistent with the 



FWZI > 6000 kms 1 measured by Filippenko (1987) . The equivalent widths of the broad component, 



relative to the modelled power-law (i.e. quasar) component of the continuum in the regions C, E, F and 
G are also shown in Ta ble [|. The fact that these are lower than the values found in typical quasars (~ 250 
A, Espey et al. (1989)) suggests that what we have modelled as the "power law" is not purely scattered 



quasar light, but may include a contribution from a non-quasar source, such as stellar populations with 
ages younger than those we have assumed in our fit. 

Several features of these emission line data are worth highlighting: 

• Apart from region A which shows a significantly lower ionization state, there is a remarkable 
consistency of lin e ratios and ioniza tion between every spatial aperture. This lack of variation 



was first noted by |Filippenko (1987)| who discussed the integrated spectra of the main NW and SE 



components and noted their general similarity. 

A high ionization state is indicated by unusually strong HeIIA4686, [NeV]A3426 in comparison with 
the nuclear NLR of most other powerful radio galaxies (see Section 11, Figure ||), and by the 



presence of the higher-ionization lines, [FeVII]A6087 and [FeX]A6375, in regions E, F and G. 
Despite the high ionization state, the temperature-sensitive [OIIIJA4363/A5007 ratio is smaller than 



measured in the nuclear and ext ra-nuclear NLR of most other powerful radio galaxies ( Tadhunter 



et al. 1989b; Binette et al. 1996). 



3.5 Physical conditions and ionization parameter 

Assuming that the AGN is the dominant photoionizing source, and that the line-emitting clouds are 
optically thick to the photoionizing radiation, values for the ionization parameter, U, (defined as the 
ratio of the number density of photoionizing photons to the number density of photoionized electrons) the 
electron density, n e> and the electron temperature, T e , may be derived from the [OII]A3727/[OIII]A5007, 
[SII]A6716/A6731 and [OIH]AA(5007 + 4959)/4363 line ratios, respectively. 

The temperat ure and density were calculated using the Lick 5-level atom calculator, fivel, flDe| 

Robertis et al. 1987| ) which, given a nominal density and value of the temperature diagnostic line ratio, 



can be used to give the temperature of the line-emitting region, and vice versa. The results of this process, 
using the de- reddened line r atios, are shown in Table |(| 

Fenston et al. (lWTfgive the following relationship between [OII]A3727, [OIII]A5007, and U: 

S§§ 

\ This relationship was calculated using a set of models in which optically thick clouds, with a constant 
gas density of 100 cm~ 3 and solar abundances, were photoionized by a central source. [OII]/[OIII] was 
chosen because it is known to be insensitive to the detailed continuum shape, while retaining sensitivity 
to U over a wide range of values. 

Values of U were calculated for each region together with the density and projected distance — 
measured from the mean continuum centroid of the SE nucleus — and these results are shown in Table 

1 

Again the results presented in Table |6j show a remarkable consistency in physical conditions and 
ionization between the various regions. 
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4 Testing the ionization models for 3C 321 



The results presented so far emphasize the fact that it is reasonable to assume that the ionization of 
the emission line gas in 3C 321 is dominated by AGN photoionization, with no clear evidence for shock- 
ionization. The justifications for this are: 

• Kinematics. The kinematics of the line emission gas are consistent with purely gravitational 
motion, and imply that the gas is relatively undisturbed. 



Phys ical conditions. Shock m odels which include a combination of cooling and precursor gas 



(e.g., [Dopita fc Sutherland 1996| ), predict relatively high electron temperatures (T > 15,000K), 



whereas the temperatures measured in several apertures across 3C 321 are much lower than this. 

Although there is no evidence for shock ionization, shocks may nonetheless influence the physical con- 
ditions in 3C 321. Such influence may be indicated by the consistently high densities measured across 
the central regions of 3C 321 (see Table ^). Warm emission line clouds in pressure equilibrium with the 
hot X-ray haloes of isolated elliptical galaxies and groups of galaxies would have electron densities in the 
range n e < 10 cm -3 at a similar radius ( Forman et al. 1985| ; Ponman et al. 1994)]. The larger densities 



measured in 3C 321 suggest a higher confining pressure, perhaps associated with the cocoon of the radio 
source. Indeed, the cocoons of radio sources are predicted to be substantially over-pressured with respect 
to the ambient ISM. 

Alternatively, the warm photoionized clouds, initially over-pressured with respect to their sur- 
roundings due the large temperature increase caused by photoionization, may not yet have expanded to 
reach pressure equilibrium with the hot ISM. 

4.1 Photoionization models 

Initially, simple photoionization models comprising a power-law continuum, f v oc v a " , photoionizing an 
optically thick, solar-abundance medium, were used to predict values of various emission line flux ratios, 



using the mappings photoionization code of |Binette et al. (1997)| . A selection of diagnostic diagrams is 
shown in Figure [j| Two power-law sequences are plotted, one with a v fixed at —1.5 and 10~ 4 < U < 0.1, 
and the second with U fixed at 0.01 with —2 < a„ < —1. Note that an ot v — —1.5 power-law continuum 



has been shown to give good fits to the s tronger emission lines in low redshift radio galaxies ( [Ferland 



fc Qsterbrock 1986| ; [Robinson et al. 1987| ). Also plotted for comparison are the emission line ratios for 
the nuclear NLR of various powerful radio galaxies from the literature. While these simple power-law 
photoionization models are in reasonable agreement with the line ratios measured in 3C 321 on some of 
the diagnostic diagrams, there are several problems. These include the following. 

• There is a general lack of consistency between the various diagnostic diagrams, in the sense that, 
even though the points may fall close to the power-law U sequence on some of the diagrams, the 
points move to different positions along the U sequence between the different diagrams. 

[NII]A6583/[OII]AA3727 is under-predicted by the power-law models (Figure ||c). This problem has 



been noted before in the case of Cygnus A ( fTadhunter et al. 1994| ) where it was ascribed to an 
N/O abundance ratio which is enhanced relative to the solar value. On the basis of Figure ||c it 
seems safe to conclude that the N/O abundance ratio is also enhanced in 3C321 and several other 
powerful radio galaxies, given that the [NII]A6583/[OII]A3727 ratio is relatively insensitive to the 
ionizing continuum shape, ionization parameter and cloud column density. 

HeIIA4686/H/3 is under-predicted by the power-law sequence (see Figures §a, b & e). This is a 
particular problem, since the ratio is not sensitive to abundances, physical conditions or ionization 
parameter. 

[OI]A6300/[OIII]A5007 and [NeIII]A3869/[NeV]A3426 are also not well-fitted by the power-law se- 
quence (see Figures |S|e, f & g), falling at or off the extreme high ionization end of the U sequence. 



2n 



The X-ray lumin osity of 3C 321 (< 10 43 erg s 1 ) is in the range typical of isolated ellipticals or groups of galaxies 



(Fabbiano et al. 198^) 
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• The locus of points formed by 3C321 and the other radio galaxies in HeIIA4686/H/3 vs 

[NeIII]A3869/[NeV]A3426 diagram (Figure ||e) is almost orthogonal to the power-law U-sequence, 
with 3C321 at the extreme high ionization end of the sequence. 

Next, in an attempt to get a better fit to the high ionization lines, photoionization models were con- 
structed for optically thick clouds photoionized by a black body spectrum with a range of characteristic 
temperatures and ionization parameters. Again two different sets of models are shown in Figure |. For 
one set, a black-body temperature of Tb = 160, 000 K was assumed, and the ionization parameter was 
varied over the range 10~ 4 < U < 0.1, while for the other set U was fixed at 10 -2 and the temperature 
varied in the range 60, 000 K < T b < 200, 000 K. Although the T b = 160, 000 K sequence gives better fits 
to the observed HeIIA4686/H/3 and [OI]A6300/[OIII]A5007 ratios, the [NeIII]A3869/[NeV]A3426 ratio is 
still a problem and, as with the power law models, there is a general lack of consistency between diagrams 
in the position of the points along the U sequence. 

We conclude that the optically thick models cannot fit the ratios of all the emission lines in a 
consistent fashion between the different diagnostic diagrams, regardless of the shape of the ionizing 
continuum. In the light of this conclusion, we consider a third set of the models which include a mixture 



of optically thick and optically thin clouds. Developed by [Binette et al. (1996)| in order to explain the 
discrepant high ionization line ratios and high electron temperatures observed in many active galaxies, 
these models are shown in Figure ^| as sequences in A M /j, the ratio of the solid angle covered by the 
optically thin component relative to that covered by the optically thick co mponent. Covering the range 
0.04 < A M /i < 16, they were generated using the tabulated values from |Binette et al. (1996)| for an 
illuminating power law with a v = — 1.3. As can be seen from Figure ||, the mixed medium models 
provide a remarkably good fit to the observed line ratios on most of the diagnostic diagrams, with 
a consistent position of the points along the A M /j sequence in all diagrams. Particularly notable is 
the fact that the A M /j sequences provide an excellent fit to the locus of radio galaxy points on the 
[NeIII]A3869/[NeV]A3426 vs HeII/H/3 plot (Figure|e), which is problematic for the optically thick models. 
In terms of the diagnostic diagrams, the only significant problem with the mixed medium models is that 
they over-predict the temperature-sensitive [OIII]A4363/A5007 ratio. However, better agreement with 



this line ratio can be obtained with super-solar abundances (see Figs. 12 & 13 of Binette et al. (1996)). 

A general problem for all photoionization models is the observed uniformity of the line ratios across 
the emission line nebula. In the case of the optically thick models, such uniformity would imply a small 
range in ionization parameter which in turn would require a 1/r 2 dependence in the density, where r is 
the distance from the illuminating source. However, such a density dependence seems unlikely, given that 
the density, where it can be measured, remains constant over a range of projected radii, and also given 
that one of the emission line regions (C) is associated with the secondary nucleus. 

The uniformity may be even more of a problem for the mixed medium models, since it requires 
not only a small variation in A M /j between the different regions, but also consistent column depths for 
the optically thick and thin components and consistent properties for the photoionizing continuum. In 
the context of the mixed medium models, such uniformity may imply some as yet unrecognized physical 
mechanism which maintains a "natural" A M / T ratio over a range of radii. 

4.2 Photoionizing photon luminosity, Q 

Another way of testing the photoionization models is to calculate the value of Q, the flux of ionizing 
photons per unit solid angle and check for consistency with other measures of the AGN brightness. Q 



may be calculated from the various physical quantities estimated in Section according to th e formula 
Q = cn e Ud 2 , where d is the (projected) distance of the line-emitting region from the source ( Robinson 



at al. 1987[ ) 



As in the previous section, we begin by considering results derived from photoionization models 
for optically thick clouds. Table |^ shows values of Q calculated for each region, under the assumption 
that the projected distance measured from our spectra correspond to the real AGN-cloud distances. The 
results reveal two problems associated with the values of Q calculated for these optically thick models: 

1. Absolute value. Q is much higher than in other low-redshift radio-loud objects - up to 2 orders 
of magnitude greater than the value measured for Cygnus A (Q ~ 10 54 photon s^sr^ 1 (|Tadhunter| 



g 



et al. 1994)), 
3C 273 (Q 



and ~ 3 times greater than the most powerful radio-loud quasar in the local uni verse 
3 x 10 55 photon s _1 sr _1 (based upon a linear interpolation across the EUV gap (Elvis 



et al. 1994|).) This is surprising in the context of the correlation between optical and radio properties 



2. 



found by Rawlings fc Saunders (1991)| since 3C 321 is not an unusually luminous radio source 
PnsMHz is roughly 200 times less than for Cygnus A (Spinrad et al. 1985). Note that the values 
of Q in Table || are lower limits, since the distances used in the calculation are projected rather 
than true distances. Region C gives the most reliable estimate of Q because it is the region with 
the least fractional uncertainty in distance from the putative illuminating source; it also gives the 
largest value of Q. 

Predicted far-infrared flux. Using the value of Q it is possible to predict the far-infrared 
luminosity, assuming that all the photoionizing photons incident on the dusty torus are reprocessed 
into the far infrared. If we assume a cone opening half-angle of 50°, the solid angle subtended by 
the torus is 47rcos50 = 8.1 sr. Assuming an isotropic illuminating continuum with the form of 
a 160,000 K black body, as indicated by the photoionization models above, with a mean photon 
energy of e — 42 eV, then the total bolometric lumin osity expected in the fa r infrared from dust- 
reprocessing is Lir = Q x e x 4ir cos 50 = 5 x 10 39 W. Golombek et al. (1988) measure the observed 
25 fim and 60 /im fluxes of 3C 321 as 353 and 1067 mJy respectively. Assuming that a power-law 
obtains between 25 and 60 //m, the integrated rest-frame luminosity between these wavelengths is 
^25-60 = 2 x 10 38 W - a factor of ~ 25 less than that expected from the measured value of Q. 



By contrast, if the mixed medium models are correct, Binette et al. (1996) estimate that the 'true' 
ionization parameter in the optically thick clouds, from which almost all of the [SII] density diagnostic 
emission originates, is U ~ 5 x 10~ 4 . The last column of Table || shows Q for the different regions 
calculated using the value of U appropriate for a mixed medium model. Repeating the above calculation 
for the expected re-processed far infrared with the values of Q calculated for a mixed medium gives 
Lin ~ 2.5 x 10 38 , which is close to the value inferred from the observations of Golombek et al. (1988) 



Although there is some uncertainty about the exact value of U for the mixed medium model, the estimates 
of Q based on this model are clearly more realistic than those derived for the optically thick case, although 
region C, with the most reliable estim ate of Q, still gives a value ~ 5 times greater than that estimated 
for Cygnus A ( Tadhunter et al. 1994 ). This is perhaps not too surprising since it has been shown tha t 
there is a large degree of scatter in quasar power for a given extended radio power ( Tadhunter et al. 1998|) . 
Also, it has been suggested that Cygnus A has enhanced radio emission in re lation to its emission lin e 
luminosity, probably as a result of its occupying a dense cluster environment ( Barthel fc Arnaud 199G| ) . 
By comparison, 3C 321 is an isolated galaxy. 

The more realistic estimates of Q for the mixed medium models further strengthen the case of 
these models in comparison with the models which assume optically thick clouds. 



4.3 The 'Q effect' 

Whether optically thick or mixed medium models are assumed for the calculation of Q, the resulting values 
show a prominent dip across the object, which coincides with the assumed position of the continuum 
source. As Q is the ionizing luminosity 'seen' by the clouds, it should be constant for all positions, since 
all clouds should 'see' the central source. This 'Q effect' has been reported before in the case of Seyfert 



galaxies ( Metz 1997 ; Metz et al. 1997 ). These authors noted that the size of the region across which the 
dip in Q was most acute (corresponding with our region F) was on the same scale as the seeing disc. This 
led them to infer the presence of an unresolved dense core, or inner narrow line region (INLR), situated 
between the conventional broad and narrow line regions. The HST [OIII] image (Figure ||) provides some 
evidence that this may also be the case in 3C 321, with the brightest pixels bordering the dust-lane in 
the SE nucleus. However, the consistency of our emission line data across the object show no evidence 
for such a dense inner region, and the highest electron densities are found in the regions furthest from 
the centre. 

A second, more plausible effect is suggested by our data - owing to the consistency of the emission 
lines across the nebula, most of the variation in Q is due to the increase in measured distance from the 
continuum centre (Table ^). These distances are clearly projected distances, and therefore represent lower 
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limits to the true values. There may therefore be a projection effect at work, whereby the line-emitting 
clouds in regions E & F are in fact some distance away fro m the AGN a l ong our line of sight. With 
respect to the unified schemes of radio galaxies and quasars ([Barthel 198$ [Lawrence 199l| ), this would 
imply that we must be viewing the central regions of 3C 321 close to the edges of the photoionized cones 
(any further towards the centres of the cones would result in a direct view of the quasar nucleus). In 
this context, it interesting to note that infrared observations by Tadhunter et al. (1999) show that the 
emission cones in Cygnus A are hollowed out, with much of the line emission generated at the edges of 
the cones. This suggests that the the projection effects indicated by the dip in Q may plausibly occur 
without exposing the hidden quasar nucleus. The current bes t estimates of cone opening angles, derived 
from examination of the statistics of AGN ( Lawrence 1991 ), therefore imply inclinations of the radio 
axis with respect to the plane of the sky in the range ~ 30 — 45° if the edge of the cone in 3C 321 is 
close to the line of sight. In this case we would expect the core radio flu x to be boosted with resp ect to 
that of the lobes, as is the case for broad line radio galaxies (BLRGs). [Vlorganti et al. (1997a) derive 
median values of of Ri.zghz = 0.0269 and 0.00295 for broad- and narrow-line radio galaxies respectively, 

Where R 2 .3GH z = P2.3G Hz, core/ '( P 2.3G Hz, total ~ Pl.ZGH z.core) ■ The fact that R 5 GHz = 0.0238 for 3C 321 

( Nilsson 1998 ) suggests that these projection effects may indeed be important, as this is clearly closer to 
the median BLRG value than that of NLRGs. This is also supported by the one-sided core-jet structure 



to the NW, visible in the maps of Baum et al. 



In summary, it is likely that projection effects can account for the dip in Q at the SE nucleus and 
that the true value of Q is close to that measured in region C, the furthest spatial bin from the nucleus 
in which we can measure Q. 



5 Conclusion 

We have presented high quality, long-slit spectral data for 3C 321. The kinematics, emission-line ratios 
and physical properties show an unusually small variation across the spatial extent of the object, and 
provide good evidence that 3C 321 is mainly photoionized, with little or no contribution from shock 
ionization. 3C 321 is therefore an excellent source to use to test the photoionization models in depth. 

We have carried out accurate continuum modelling and subtraction to allow the measurement of 
faint emission lines so as to better constrain the photoionization models, and allow calculation of physical 
conditions across the nebula. The continuum modelling requires the presence of a significant, spatially 
extended, young stellar population to give good fits to the data. 

Comparison of emission-line flux ratios with the photoionization models shows that the mixed 
medium models - involving a combination of clouds which are optically thick and thin to the ionizing 
radiation - provide an excellent fit to the observed line ratios, especially if super-solar abundances are 
assumed. Conventional models, which assume a single population of optically thick clouds, do not provide 
a good fit to the observed line ratios, and require differing values for parameters such as U to reproduce 
the line ratios in different diagnostic diagrams. 

Our results also show that caution is required when using the emission line properties to calculate 
the ionizing luminosity of the illuminating AGN (Q). Under the assumption that the line emitting clouds 
are optically thick, we find that the estimated value of Q is unrealistically large, given the radio and far- 
infrared properties of the source. However, the lower value of ionization parameter appropriate for the 
mixed medium model leads to an estimate of Q which is more consistent with the estimated far-infrared 
luminosity of 3C 321. 

The unusual degree of consistency of the emission line ratios across the nebula presents difficulties 
for any ionization model due to the need for consistent physical conditions and ionization over large range 
of distances from the illuminating nucleus. This is a particular problem for the mixed medium model 
which provides the best fit to the line ratios, since the uniformity requires the ratio of covering factors of 
the optically thick and thin components to remain remarkably constant across the nebula. 
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Region 


15 Gyr flux T 


Power law flux T 


a 


■ 


Goodness-of-fit* 


A 
B 
C 
D 
E 
F 
G 


67+j 
80±li 

69 ±7 
71 +io 

78+.? 


271* 
31±g 

23+5 
27±1 
30+| 
20t^ 


u -°-1.0 

-1.3 ±0.6 

-0.9+25 

-i o+0.4 
J -°-0.6 
_i 1+0.6 
1 - i -0.8 

o i+°- 4 

u,1 -0.6 

_n 9+ 6 

u - y -0.8 


10.98 
2.45 
0.77 
4.09 
1.44 
0.55 
1.96 


0.000 
0.000 
0.820 
0.000 
0.051 
0.982 
0.001 



Table 1: 2-component fits to the continuum in the each of the regions A-G. Errors quoted denote the 
range of values which give % 2 within the range Xmin < X 2 < Xmin + A(x 2 ), where A\ 2 is the interval in 
X 2 corresponding to the 68% confidence level (Press et al. 1992). The value of Xmin red measured for A is 
almost certainly an overestimate - the error on each data point is much greater than the flux calibration 
error due to noise in the continuum. 

'As a percentage of the measured continuum after subtraction of the nebular component in the rest 
wavelength region 5046.0- 5075.0 A. 

^Defined as the formal statistical probability that this value of Xmin red nas not Deen arrived at by chance. 



Region 


15 Gyr flux' 


0.1 Gyr flux T 


Power law flux T 


a 


%min,red 


Goodness-of-fit 


A 
B 
C 
D 
E 
F 
G 


62-5™ 
57.5 ± 12.5 

79 r + 12.5 

' z -°-17.5 

52.5 ±5.0 
55 0+ 17 - 5 
60.0+f 5 5 
67.5+1° 


Q+b 

20.0± 7 2 5 5 
i n n+ 12 - 5 

1U.U_ 10 Q 

40.0+J;° 
25.0 ± 12.5 
15 0+ 75 

io - u -15.0 
9(^+12.5 
zu -7.5 


27.5 ± 10.0 

99 r+10.0 
. o 27 5 

20.0 ± 12.5 

o o+ 7 - 5 
u - u -o.o 

17 5+ 7 - 5 
25.0+ 7 2 5 5 
1 2 5+ 7 ' 5 


-o.9±iS 

n/a 

-o.i±S:8 

0.9t°;« 


11.32 
1.24 
0.64 
0.74 
0.38 
0.44 
1.12 


0.000 
0.139 
0.921 
0.818 
0.999 
0.995 
0.253 



Table 2: As Table [j], but for three components. The power law flux in region D makes no contribution 
to the best fit model and so all values of a give an equally good fit. 

^As a percentage of the measured continuum in the rest wavelength bin 5046.0 — 5075.0 A. 



Region 


Mass in 15 Gyr component 
/10 10 M Q 


Mass in 0.1 Gyr component 
/1O 8 M 


Percentage of young 
stars by mass 


A 


0.23 






B 


0.75 


0.3 


0.4 


C 


4.19 


0.6 


0.1 


D 


3.38 


2.6 


0.8 


E 


3.79 


1.7 


0.5 


F 


4.98 


1.3 


0.3 


G 


2.52 


0.8 


0.3 


Total 


19.9 


7.2 


0.4 



Table 3: Total mass of stars in each component. These values are obtained by scaling the mass of the 
original models (1 M Q ) to give the luminosity (L /A) in the normalizing bin indicated by the model fits. 
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Probable 
ID 


A 


B 


Measured values relative to H/3 
C D E 


F 


G 


[NeV]A3346 





36 ±2 


27 ± 1 


24 ± 2 


46 ± 3 


37 ±1 


35 ±1 


[NeV]A3426 


50 ± 10 


100 ±6 


75 ± 2 


65 ± 6 


127 ± 9 


104 ±3 


97 ±3 


[OII]A3727 a 


229 ± 12 


131 ±4 


128 ± 2 


118 ± 5 


141 ± 6 


111 ± 2 


98 ±2 


[NeIII]A3869 


77 ±9 


77 ±4 


69 ±2 


55 ±7 


88 ±8 


68 ±2 


62 ±2 


[NeIII]A3964 b 


24 ±6 


32 ±3 


30 ±2 


45 ±6 


24 ± 7 


32 ±2 


32 ±2 


[SII]A4069 a 


12 ±5 


6±2 


5± 1 






4± 1 


4.3 ± 1.0 


R5 


26 ±5 


26 ±3 


14 ± 1 


18 ±3 


13 ±3 


18 ±1 


17 ±1 


H7A4340 


34 ±5 


43 ±3 


40 ±1 


51 ±5 


37 ±4 


36 ±2 


39 ±4 


[OIII]A4363 


— 


10 ±3 


7 ± 1 


5.3 ± 3.2 


9 ± 3 


8± 1 


13 ±3 


HeIIA4686 


— 


50 ±3 


51 ± 1 


47 ± 4 


54 ± 3 


45 ±1 


47 ±1 


([ArIV])A4713 





6±2 


8 ± 2 






9± 1 


8±2 


H/3A4861 


100 


100 


100 


100 


100 


100 


100 


[OIIIJA4959 


288 ± 14 


305 ±7 


288 ± 4 


301 ± 10 


366 ± 13 


356 ±6 


363 ±6 


[OIII]A5007 


860 ± 41 


910 ±21 


902 ± 11 


899 ± 30 


1094 ± 40 


1062 ± 17 


1083 ± 19 


[NI]A5199 








5 ± 1 


< 6 




6± 1 


5± 1 


[FeVII]A6087 












_ 


22 ±2 


13 ±1 


[OI]A6300 





10 ±2 


17 ±3 


— 


17±5 


18 ±1 


18 ±1 


[OI]A6364 


— 


< 4 


6± 1 






6± 1 


6± 1 


[FeX]A6375 










8±3 


7± 1 


4±2 


[NII]A6548 


54 ±3 


47 ± 1 


70 ±2 


52 ±2 


77 ±3 


82 ±2 


57 ±3 


HaA6563N 


406 ± 22 


321 ±8 


381 ±6 


283 ± 10 


398 ± 15 


483 ±8 


434 ±11 


HaA6563B 






35 ± 14 




88 ± 14 


115 ± 16 


91 ±24 


[NIIJA6583 


163 ± 10 


140 ±4 


212 ±5 


155 ±7 


231 ±9 


247 ±6 


170 ±8 


HeIA6679 




< 2 


8±3 


< 2 




7± 1 


6±2 


[SII]A6716 


88 ± 13 


46 ±2 


60 ±1 


40 ±3 


63 ±3 


75 ±1 


58 ±1 


[SII]A6731 


69 ±9 


31 ±2 


57 ±1 


40 ±3 


57 ±3 


63 ±1 


50 ±1 


[ArIII]A7136 




16 ±2 


24 ±2 


19 ±3 


39 ±7 


33 ±2 


31 ±2 


H/3 flux (10~ ib 
ergs _1 cm~ 2 A _1 ) 


0.24 ±0.01 


0.51 ±0.01 


1.63 ±0.02 


0.56 ±0.02 


0.53 ±0.02 


1.54 ±0.02 


1.23 ±0.02 


H/3 equivalent 
width (A) 


29 ±7 


15 ± 1 


11.1 ±0.2 


3.6 ±0.1 


2.6 ±0.1 


7.6 ±0.2 


13.0 ±0.4 


Broad Ha 
FWHM (kms- 1 ). 






4000 ± 1000 




3300 ± 400 


4300 ± 800 


3400 ± 800 


Broad Ha 
equiv. width 

(A). 






46 ±34 




28 ± 20 


59 ±25 


200 ± 150 


E(B-V) c 


0.21 ±0.02 





0.31 ±0.01 





0.38 ±0.04 


0.50 ±0.01 


0.41 ±0.01 



Table 4: Measured emission line fluxes as a percentage of flux in H/3. The narrow and broad components 
of Ha are denoted by N and B respectively. All lines in the wavelength ranges covered by the reduced 
spectra are shown. All quoted errors are Gaussian-fitting errors, and do not include a contribution from 
flux calibration errors (estimatd as 5%). Also shown are the equivalent widths of H/3 relative to the total 
continuum and of broad Ha relative to the power-law component of the continuum model in each region. 
"Values for blends, 
includes contribution from He 

c Weighted means deduced from Ha/H/3, H7/H/3 and H<5/H/3 and assuming Case B. 
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Probable 
ID 


A 


Reddening corrected values relative to H/3 
B C D E F 


G 


[NeV]A3346 


— 


36 ± 2 


40 ± 1 


24 ±2 


74 ±5 


70 ±2 


59 ±2 


[NeV]A3426 


64 ± 13 


100 ±6 


107 ±3 


65 ±6 


197 ±13 


185 ±6 


157 ±5 


[OII]A3727 a 


274 ± 15 


131 ±4 


166 ±3 


118 ±5 


194 ±9 


172 ±4 


141 ±3 


[NeIII]A3869 


90 ± 11 


77 ±4 


86 ±2 


55 ±7 


117 ± 10 


99 ±3 


85 ±2 


[NeIII]A3964 f) 


28 ±6 


32 ± 3 


37 ± 2 


45 ± 6 


31 ± 9 


45 ±3 


42 ±2 


[SII]A4069 a 


14 ±6 


6±2 


6± 1 






5± 1 


6± 1 


H<5 


29 ±5 


26 ±3 


17 ±1 


18 ±3 


16 ±4 


24 ±1 


22 ±1 


H7A4340 


37 ±5 


43 ±3 


45 ± 1 


51 ±5 


43 ±5 


45 ±2 


47 ±4 


[OIII]A4363 


— 


10 ±3 


8±1 


5.4 ±3.2 


10 ±3 


9± 1 


15 ±4 


HeIIA4686 


— 


50 ±3 


53 ±1 


47 ±4 


57 ±4 


48 ±1 


50 ±1 


([ArIV])A4713 


— 


6±2 


8±2 






9±2 


8±2 


H/3A4861 


100 


100 


100 


100 


100 


100 


100 


[OIII]A4959 


283 ± 14 


305 ±7 


281 ±4 


301 ± 10 


355 ± 13 


341 ± 6 


352 ±6 


[OIII]A5007 


839 ± 40 


910 ±21 


867 ± 11 


899 ± 30 


1042 ± 38 


996 ± 16 


1031 ±18 


[NIJA5199 


— 




5 ± 1 


< 6 




5 ± 1 


4± 1 


[FeVII]A6087 


— 










14 ±1 


9± 1 


[OI]A6300 


— 


10 ±2 


13 ±2 




12 ±3 


11 ± 1 


12 ±1 


[OI]A6364 


— 


< 4 


4.1 ±0.6 







4± 1 


4± 1 


[FeX]A6375 


— 








5±2 


4± 1 


2± 1 


[NIIJA6548 


43 ±3 


47 ±1 


50 ±1 


52 ±2 


51 ±2 


47 ±1 


36 ±2 


HaA6563N 


323 ± 17 


321 ±8 


270 ±4 


283 ± 10 


261 ± 10 


278 ±5 


275 ±7 


HaA6563B 






25 ±10 




58 ± 10 


66 ±9 


57 ± 15 


[NII]A6583 


130 ±8 


140 ±4 


150 ±3 


155 ±7 


151 ±6 


141 ±3 


108 ±5 


HeIA6679 




< 2 


6±2 


< 2 




6±2 


4± 1 


[SII]A6716 


69 ± 10 


46 ±2 


42 ±1 


40 ±3 


40 ±2 


42 ±1 


36 ±1 


[SII]A6731 


54 ±7 


31 ±2 


40 ± 1 


40 ±3 


36 ±2 


35 ± 1 


31 ± 1 


[ArIII]A7136 




16 ±2 


16 ± 1 


19 ±3 


23 ±4 


17 ±1 


18 ± 1 



Table 5: As Table [| above, but corrected for reddening using the values of E(B-V) calculated for each 

region. Note the consistency across the different spatial regions for any given line ratio. 

a Values for blends. 

includes contribution from He 





T 




n e 


U 


Q 


Q (U = 5 x 10" 4 ) 




(K) 


(10 20 m) 


(cm' 3 ) 


(10- 2 ) 


(phot s _1 sr _1 ) 




A 




5.4 ±0.7 


< 380 


0.56 ±0.04 






B 


12000™ 


4.0 ±0.7 


< 60 


1.20 ±0.05 






C 


noootl°° 


2.5 ±0.6 


500+™ 


0.95 ±0.02 


(8.8 ±2.4) x 10 55 


(4.7 ± 1.2) x 10 54 


D 


10000™ 


1.4 ±0.6 


550+ 240 


0.71 ±0.04 


(4.2 ± 2.4) x 10 55 


(1.5 ±0.9) x 10 54 


E 


11600™ 


0.4 ±0.3 


390+ 170 
oau_ 130 


0.98i°;°g 


(2.1 ± 1.9) x 10 54 


(1.1 ± 1.0) x 10 53 


F 


11300 ±600 


0.4 ±0.5 


260 ± 40 


1.05 ±0.03 


(1.1 ± 1.4) x 10 54 


(5.3 ±6.6) x 10 52 


G 


1330011?°° 


1.3 ±0.5 


310 ±60 


1 q+0.04 
i -°-0.03 


(2.1 ±0.8) x 10 55 


(7.7 ±3.1) x 10 53 



Table 6: Table showing reddening corrected values of temperature, distance from 'nucleus,' electron 
density, ionization parameter and photoionizing luminosity, respectively. The upper limits for the density 
in regions A & B are la limits. The uncertainties in Q are dominated by the uncertainties in the distances, 
which correspond to the width of the regions A-G. Also shown is the value of Q derived for each region 



assuming a value of U more appropriate to the mixed medium models discussed in Section 4.1 
^Assumes H = 50 kms _1 Mpc~ 1 , q Q = 0.5 
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Figure 1: (See accompanying jpeg file.) Position 
image taken using the HST. The different regions 
detailed in Section 3 were carried out are labelled 
(jVlartel et al. 1999[). 



of the slit superimposed on an [OIII]AA4959 + 5007 
for which the various measurements and calculations 
A-G. Note that there is a gap between regions B & C 



i • • • • r 




Shift from continuum centre (arcsec) 



Figure 2: Plot showing the spatial variation of [OIII]A5007 emission (solid line), together with that of 
a continuum bin centred on 5110 A (dashed line) in the rest frame. The positions of the spatial bins 
chosen for analysis are also shown. The continuum flux has been scaled by a factor of ~ 50 for ease 
of comparison. '0' is defined as the position of the continuum measured by fitting a Gaussian to the 
continuum spatial profile in a bin centred on 8849 A. 
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Figure 3: ISIS spectra of the different regions of 3C 321, together with best-fit composite 3-component 
model spectra and individual model components. Individual plot legends provide the key to the com- 
ponents displayed. The 3-component model which minimizes x 2 in each case can be seen to fit the 
continuum in the different regions extremely well. 
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Figure 3: Continued 
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Figure 4: Emission-line kinematics. Plots tracing the variation in a) line centre and b) line width across 
3C 321. The solid lines denote the 1 sigma error bars on [OIIIJA5007 measured on a pixel-by-pixel basis 
across the image. The triangles denote Ha, diamonds H/3 and squares [SIIJA6717. 
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a) 





Figure 5: Diagnostic diagrams showing a) [OIIIJA4363/5007 vs HeIIA4686/H/3, b) 
[OII]AA3727/[NeV]A3426 vs HeIIA4686/H/3, c) [NII]A6583/[OII]AA3727 vs [OII]AA3727/[OIII]A5007, 
d) [OIII]A4363/5007 vs [OII]AA3727/[OIII]A5007, e) [NeIII]A3869/[NeV]A3426 vs HeII/H/3, f) 
[OII]AA3727/[OIII]A5007 vs [OI]A6300/[OIII]A5007, g) [OIII]A5007/H/3 vs [OI]A6300/[OIII]A5007, 
h) [SII]A6724/Ha vs [OII]AA3727/[OIII]A5007, i) [OIII]A5007/H/3 vs [OII]AA3727/[OIII]A5007. Dashed 
lines represent power-law ionizing continua, / oc v a " , with triangles representing a v — —1.5, U = 1CP 4 , 
2.5 x 10~ 4 , 5 x 1(T 4 , 10~ 3 , 2.5 x 1CT 3 , 5 x 10" 3 , 0.01, 0.025, 0.05, 0.1 and crosses (x) showing models 
with U fixed at 0.01 and a v = —2.0, —1.75, —1.5, —1.0. Dotted lines are black body models, with 
squares representing sequences of models with T — 160, 000 K and the same sequence in U as before, and 
crosses showing sequences with U = 0.01 and T = (033, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0) x 10 5 K. Solid lines 
correspond to models with a«„ = —1.3 power law illuminating mixture of optically thin and thick clouds. 
The reddening corrected values for 3C 321 are shown as filled diamonds. '+' denotes measurements for 
other radio galaxies taken from iDickson (1997)1 iTadhunter (1983)1, iGrandi & Osterbrock (1978)1, ICosterol 
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Figure 5: continued 
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Figure 5: continued 
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Figure 5: continued 
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